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Summary. The airborne traffic monitoring system forms a novel technology of 
detecting vehicle motion. An optical digital camera located on an airborne platform 
produces a series of images which then are processed to recognized the fixed vehicles. 
In this way the video data are converted into the time sequence of frames containing 
the vehicle coordinates. In the present work a three-frame algorithm is developed 
to identify the succeeding vehicle positions. It is based on finding the neighboring 
points in the frame sequence characterized by minimal acceleration. To verify and 
optimize the developed algorithm a "Virtual Road" simulator was created. Finally 
available empirical data are analyzed using the created algorithm. 

Introduction 

The traditional techniques of traffic flow measurements are based on local detec- 
tors mounted at fixed places of a road network. Such detectors provide adequate 
information about traffic flow only for road fragments of a rather simple structure. 
However in cities there are a large number of complex road intersections, where ve- 
hicles moving on different branches interact strongly with one another. To measure 
appropriately traffic flow features on such "hot areas" the information supplied by 
the stationary detectors should be complemented at least with the detailed infor- 
mation about the spatial structure of traffic streams on these hot areas. Currently 
the novel technology of airborne monitoring and detecting the main characteristics 
of traffic flow is under development. The German Aerospace Center (DLR) with 
its department of Transportation Studies and department of Optical Information 
Systems has developed an airborne traffic monitoring system which was successfully 
used during the Soccer World Championship 2007 in Germany [1, 2, 3]. In particu- 
lar, these and some other DLR airborne datasets are used in the present paper. This 
system is based on an optical sensor placed on an airborne platform. The optical 
sensor produces a series of images of area under. The empirical airborne traffic data 
obtained by recognizing these images form a set of the frames with individual ve- 
hicle coordinates. To use these data and to calculate the main parameters of traffic 
flow trajectories of individual vehicle motion have to be reconstructed. In this way 
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Fig. 1. Illustration of the car identification algorithms. Within the two-frame speed 
limitation algorithm (left fragment) new objects {x n ,y n } are sought in the circle of 
radius r — v ma xT centered at the point {a;™ -1 , ln the three-frame algorithm 

bounding acceleration (right fragment) new objects {x n ,y n } are sought in the circle 
of radius r — a max r 2 centered at the point {2a;" -1 — x™~ 2 , 2j/™ _1 — yZ~ 2 }- 

it becomes possible to measure the traffic flow rate, the mean velocity of vehicles, 
and their density simultaneously. 

Vehicle Trajectory Reconstruction 

The airborne traffic data under consideration are the time series T of frames {Ft} 
with the coordinates of recognized vehicles, T = {F t = {x a ,y a }}. The problem is 
to construct a collection of trajectories V = {x(t),y(t)}t passing through the points 
{{xc^yc}} at the corresponding time moments {t — nr}, where n is integer and r 
is the time span between successive frames. Some constraints should be imposed on 
the vehicle trajectories V to make them smooth. There are several ways to do this, 
in particular, to bound the "acceleration" (the used approach) 

or to bound the "velocity" 




We note that the acceleration and velocity thresholds, a max and « m ax, are not the 
real characteristics of car motion but internal algorithm parameters whose choice is 
determined by its efficiency being highest. 

To explain the crux of the algorithm implementing these constrains let us assume 
that the preceding frames {...,F n _2, -Fn-i} have been analyzed and the vehicle 
trajectories are reconstructed at the previous time moments {. . . ,n — 2, n — 1}. 
Then a car {x^,y2} in the frame F n will be incorporated into the trajectory ui = 
U m {x™,y™}if 

yj (xZ~ 2 - 2XZ' 1 + zS) 2 + (yZ~ 2 - 2y£- 1 + y£) 2 < a ma * r 2 . (2) 
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object loss (type A) trajectory mixing (type B) 



Fig. 2. The main errors of the car identification algorithm. 



Figure 1 illustrates the given three- frame algorithm based on (la) as well as the 
relative two-frame algorithm based on (lb) mentioned for comparison only. 

Let us discuss the possible errors of the car identification algorithm caused by 
errors in the car coordinate measurements. There are two main types of the latter 
ones; the errors of individual vehicle positions, £i, and the errors in the frame ref- 
erence to GPS, £2. According to the pilot airborne monitoring system (Institute for 
Traffic Research, DLR, Berlin) [1, 2] the error values can be estimated as £1 ~ 0.5 m 
and £2 ~ 1-2 m. The caused errors of the car identification algorithm, i.e. the ob- 
ject loss (type A) and the trajectory mixing (type B) illustrated in Fig. 2 will be 
quantified in units of 

^/ m The number of lost(mixed) objects -^qqo/ 
Total number of objects 

To analyze the feasibilities of the identification algorithm, first, the following 
kinematic simulator has been used (Fig. 3, upper fragment). The motion of cars is 
specified by the shown equations where the parameters xoi, yoi, b X i etc. are ran- 
dom and some addition constraints are imposed to prevent the collisions. The noise 
components £ x (y)i = (£1 + &) x (y)i imitate the errors in the measurements of the 
car coordinates Xi and tji at the time moments t = nr. The system parameters 
and the number of cars are specified so that the virtual car ensemble imitate the 
real characteristics of traffic flow, at least semiquantitatively. The lower fragment in 
Fig. 3 exhibits the net value of the identification errors, A + B, vs. the acceleration 
threshold a max . As could be expected, the present results demonstrate the fact that 
the car identification algorithm attains its maximal efficiency when 

_ 4(£im + 6a/) ,„s 

Umax — 2 ' \ / 

and, in particular, for £im = 0.5 m and £2 = 1.5 m the optimal acceleration threshold 
is about a max ~ 50 m/s 2 . 

Second, the developed algorithm has been verified using the available airborne 
traffic data collected during Soccer Championship 2006 in Stuttgart within the DLR 
project "Soccer" [3]. These data undergone manual processing, so initially we had a 
collection of vehicle trajectories regarded as real keeping in mind the human abilities. 
The vehicle identification algorithm has given the result depicted in Fig. 4. For 
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Fig. 3. Airborne data simulator, "Virtual Road", used to verify the vehicle identi- 
fication algorithm (upper fragment) and the corresponding results of the car iden- 
tification routine (lower fragment). The subscript M at the noise components £;m 
stands for their amplitudes. 
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Fig. 4. The airborne data collected during the Soccer Championship 2006 in 
Stuttgart. Typical example of the analyzed video data (left fragment) and the results 
of car identification routine vs. the acceleration threshold (right fragment). 



these data, again, the optimal value of the acceleration threshold meets relation (3). 
Figure 4 presents also the efficiency of the car identification algorithm for different 
values of the velocity threshold « ma x. As seen the imposition of the additional speed 
limitation constraint reduces the algorithm efficiency. In other words, constraints (1) 
interfere with each other and, thus, should be used separately. As should be expected 
the measurement errors and the optimal acceleration threshold are again related by 
expression (3). 
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Fig. 5. Time-mean speed (green) and space-mean speed (red) reconstructed using 
the smoothed vehicle trajectories (left fragment) and the same data obtained using 
the neighboring points directly. The airborne data of pilot flights in Berlin, 2004, 
the LUMOS project. 
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Fig. 6. Fundamental diagram (right fragment) and an example of the analyzed 
video data (left fragment). The NGSIM "First Prototype" dataset [4]. 



Mean Velocity of Traffic Flow and Fundamental Diagram 

Because of the considerable errors in the vehicle coordinates the found successive 
positions of a car cannot be used directly to calculate its velocity. However, the 
reconstructed trajectories enable one to overcome this problem by fitting a rather 
smooth trajectory to the found data. In this way we have analyzed the empirical 
data collected during pilot flights of an airplane in Berlin, 2004 (DLR, LUMOS 
project). The mean velocity of traffic flow was calculated, first, by smoothing the 
reconstructed vehicle trajectories via the Savitzky-Golay filter and, then, averaging 
the velocities over the observed car ensemble. The results are illustrated in Fig. 5. As 
seen the time pattern of traffic flow speed obtained by smoothing the reconstructed 
vehicle trajectories gives much more adequate description of traffic flow. 
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Using the developed technique we have analyzed also the NGSIM "First Proto- 
type" dataset consisting of vehicle trajectories on a half-mile section of Interstate 
80 in Emeryville, California, for one half of hour [4]. The system portrait on the 
phase plane "traffic flow rate - car density" was drawn using the mean car velocity 
calculated via the developed technique and car density obtained by just counting 
the detected objects (Fig. 6). No widely scatted states are visible, instead, branching 
is fixed. This result again poses the question about the necessity to single out the 
vehicle flow discreteness in analyzing the states of traffic flow [5] . It should be noted 
that the developed approach depresses this discreteness effect. 

Conclusion 

We have developed the three-frame-algorithm of reconstructing vehicle trajectories 
from airborne traffic data, which evaluates the proximity of objects in the neighbor- 
ing frames in terms of acceleration. The acceleration threshold a max is determined 
by the car measurement accuracy. Using the "Virtual Road" simulator and the 
dataset collected within the project "Soccer" (DLR, Stuttgart, 2006) the value of 
Umax ~ 40-50 m/s 2 is found and the possibility of decreasing the identification error 
down to 5-30 % is shown. 

Using one of the NGSIM datasets it has been demonstrated that for the observed 
congested traffic the fundamental diagram splits into two branches when the car 
density exceeds some value rather than exhibits the widely scatted states. This result 
again poses a question whether the discreteness of traffic flow could be responsible 
for the appearance of widely scatted states on the fundamental diagrams. Naturally, 
the unambiguous answer requires special and detailed investigation. 
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